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Purpose: Grating-based X-ray phase- contrast imaging is a promising modality increasing the 
soft tissue contrast in medical imaging. In this work, the signature of lung tissue in X-ray grating- 
based physe-contrast imaging is investigated. 

Methods: We used a Talbot-Lau interferometer for our investigations of two C57BL/6 mice. Both 
underwent projection imaging and computed tomography. 

Results: The results show that the three images obtained by X-ray phase-contrast imaging show 
complementary anatomical structures. Especially the dark field image allows a more-exact determi- 
nation of the position of the lung in the chest cavity. 

Conclusion: Due to its sensitivity to granular structures, the dark field image may be used for the 
diagnosis of lung diseases in earlier stages or without a CT scan. Furthermore, X-ray phase- contrast 
imaging may also have great potential in the application of animal laboratory sciences to reduce 
the number of required animals used in long-term translational, toxicity, and regenerative medicine 
studies. 
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FIG. 1: Sketch of the used Talbot-Lau Interferometer with X-ray source S, source grating GO, a sample O, the phase 
grating Gl, the analyser grating G2, and a pixelated detector D. I(x) illustrates the idealized, monoenergetic Talbot 

intensity pattern, produced by Gl and distorted by the object. 



I. INTRODUCTION 



Since the discovery of X-ray radiation by W.C. Rontgen in 1895, its attenuation properties have been used in many 
fields such as material science and medicine. On the other hand, as an electromagnetic wave, X-rays are not only 
absorbed when traversing matter, but also shifted in phase. To detect this phase-shift, many techniques mainly for 
the use at synchrotron facilities have been developed [IH3]. Another approach which was published the last few years 
is the usage of a grating interferometer [IHZ]. Furthermore, it has already been shown that it is possible to employ 
this technique at low-brilliance, medical X-ray sources |8J to improve the soft-tissue contrast. 

One possible medical application exploiting this ability of X-rays is lung imaging. Here, attenuation-based X-ray 
imaging is the most-often-used modality to obtain morphological information. It has been shown by Kitchen et al. 
[HI HDL that lung tissue seems to be the ideal organ for phase-contrast imaging, due to its small-sized alveoli. 

Kitchen et al. carried out their measurements at a synchrotron facility using a high flux monoenergetic X-ray beam 
and a high-resolution X-ray detector, not applicable for medical routine. In this work, in contrast, we report on 
our investigations of in situ lung tissue imaging in a grating-based X-ray phase-contrast setup using a conventional, 
medical X-ray tube and a commercial X-ray flat-panel detector as it is widely used in a medical environment. 



II. MATERIALS AND METHODS 



A. Basics of Talbot-Lau interferometry 



X-ray phase-contrast imaging can be performed with low-brilliance medical X-ray sources using Talbot-Lau inter- 
ferometry. This technique was proposed by Pfeiffer et al. [8 J . The basic idea of such a setup is shown in Figure [I] A 
medical X-ray source produces non-coherent, polychromatic X-rays. Downstream the direction of wave propagation a 
source grating (GO) ensures the spatial coherence by introducing many virtual slit sources. The wavefronts originating 
from these slit sources impinge on the object, where they become deformed depending on the object's material prop- 
erties. Further towards the detection plane a phase grating (Gl) imprints a periodical phase shift on the wavefronts 
emitted from GO. For a coherent, monoenergetic plane wave and a 7r-shifting phase grating Gl a rectangular intensity 
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FIG. 2: Example intensity oscillations I(p) detected in a detector pixel as a function of the lateral position of 
grating G2. The solid line shows the reference oscillation without the object while the dashed line shows the 
resulting oscillation with object. / are the mean values while A are the amplitudes of the oscillations. A(/> is the 
phase difference between the oscillations and is the measurement parameter for the differential-phase image. 



modulation occurs in the so-called constructive fractional Talbot distances [HJ [E] behind Gl: 
For a spherical wave this distances rescale due to the magnification of the setup to 

«=^. w 

where g\ is the grating constant of the phase grating, A the wavelength of the X-rays, / the GO to Gl distance, and n 
an integer number. 

The resulting regular intensity pattern I(x), which is in general too small to be measured directly, can be sampled 
by a measurement of the intensity for a number of grating positions p of the analyser grating (G2) called "phase steps". 
This technique was adapted by Weitkamp et at. [6 J from visible light interferometry p~3| [T4] . When this procedure 
is repeated once with and once without the sample, two intensity oscillations I(p) are obtained in each pixel of the 
pixelated detector (D). 

Figure [2] illustrates the intensity oscillations obtained in one pixel for an ideal, monoenergetic case. The subscript 
"ref" refers to the reference oscillation without object, while "obj" denotes values with the object present in the beam. 
From these two curves the three image quantities of grating-based X-ray phase-contrast imaging can be obtained: the 
absorption image as the ratio of the mean intensities J Q bj / I re f of the sampled oscillations, the differential-phase image 
Acf) as the phase difference of the two intensity oscillations, and the dark field image as the ratio of the visibilities 
V b} I Kef? with V being defined as 

t t -^max -^min A . * 

- J ZTT - ~ y W 

- 1 max r -^min 1 

of the two oscillations, where 7 max is the maximum and I m [ n the minimum detected in a pixel. This can also be 
approximated as the fraction of the oscillation's amplitude A and its mean intensity /. In the non-ideal case, where 
GO has a not-infinitely-small slit width and a photon spectrum is used, the triangular oscillations shown in Figure [2] 
become sinusoidal due to the lack of spatial coherence generated by the source grating and the fact that the phase 
shift of Gl matches tt only for one energy - the design energy - of the polychromatic X-ray spectrum. 



B. Phase retrieval 



From the two intensity oscillations I{p) recorded as described in the previous section, the three image quantities can 
be obtained in several ways. As the oscillations are sine functions, a fast Fourier transform (FFT) can be performed 
for each pixel for the reference- and the object oscillation. Alternatively, a least-squares algorithm fitting a sine 
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TABLE I: Properties of the used gratings as stated on the delivery notes. 



source 


grating GO phase 


grating Gl analyser grating G2 


grating constant (/im) 23.95 


4.37 


2.40 


thickness (/xm) 150 


8.7 


100 


duty cycle 0.50 


0.50 


0.50 


material gold 


nickel 


gold 



function to the data can be used. This technique is used for all images shown in this work. The least-squares method 
has some advantages compared to the FFT method. For example, the period of the sine function can be added as 
an additional fitting parameter to correct systematic sampling errors or grating variations due to thermal expansion. 
Furthermore, a change in the fit model function can correct for other systematic errors [15J. 

Using one of the methods mentioned above, the parameters defining a sine function the offset /, the amplitude A 
and the phase <p are determined from the measured data. These are used to calculate the three image quantities in 

each pixel: the absorption fi = =i, the differential phase A(j) = </> re f — o bj and the dark field D = y^j- 

C. Setup parameters 

All measurements presented were carried out at our benchtop phase-contrast imaging setup. It consists of a Siemens 
MEGALIX X-ray tube driven at 60 kV accelerator voltage. The anode material of the tube is tungsten and only the 
intrinsic filtering was applied. For photon detection, the Varian PaxScan 2520D flat panel detector with Csl as 
scintillation material and a pixel size of 127xl27/im 2 was used. Before the measurements were carried out, an offset 
and gain calibration of the detector was performed as recommended by the manufacturer. 

The gratings of our setup were produced by Karlsruhe Institute of Technology (KIT) employing the LIGA 
method [16J. The parameters of the gratings are summarized in Table |TJ For the calculation of the fractional 
Talbot and grating distances a design energy of 25keV was assumed. Except the source grating GO, which had a 
quadratic field of view (FOV) of 5x5 cm 2 , all gratings were of rectangular shape with an effective area of 2x6 cm 2 . 

The measurements were carried out in the second constructive fractional Talbot distance (n = 2 in equation [T]). 
Taking the source- and analyser- grating properties into account, the setup distances were calculated using equation [l] 
and the theorem on intersecting lines resulting in 161.2 cm for the GO-to-Gl distance and 15.9 cm for the Gl-to-G2 
distance. 

The focal spot was located 15.5 cm in front of the source grating GO, due to the size of the X-ray housing. Further, 
the detector was placed 5 cm behind the analyser grating, due to the dimensions of the grating mounting. The object 
was placed 10 cm in front of the phase grating Gl. 



D. Animals 

Two C57BL/6 adult male mice were randomly selected from carbondioxide-killed ex-breeding stock at the Franz- 
Penzold-Zentrum animal facility for investigation. Both animals were press fit in a 50 ml conical polyethylene tube to 
reduce any movement while the measurements were carried out. No further preparation procedures with the samples 
were done. This study complies to the institution's animal welfare and standard procedures. 



III. RESULTS 

In situ lung imaging of two healthy C57BL/6 adult mice was performed to obtain information on the signals of 
phase-contrast imaging. Figure [3] shows the three images obtained of the lungs of a mouse in posterior-anterior (p. a.) 
orientation, using eight phase steps with Ax = 0.3 um and an acquisition time of 3.3 s each at a tube current of 30 mA. 
This corresponds to an air kerma of 1.58 mGy measured at the object position with a calibrated IBA Dosimax plus 
A HV dosimeter with the solid-state detector unit RQX 70kV. 

While the lung tissue is only weakly visible in the absorption image in Figure [3] the bony structure of the mouse 
thorax can clearly be seen. In the differential-phase image the fur of the animal is clearly visible. Both lungs as well 
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FIG. 3: The three X-ray phase-contrast images of the lung region of a mouse in posterior-anterior (p. a.) orientation. 
Top: absorption image, middle: differential phase image, bottom: dark field image. 



as the spinal cord of the animal are shown. The speckle noise in the phase image described by Kitchen et at. |9J can 
be observed in the expected lung area. Dark field imaging highlights clearly the fur-air interface as well as both lungs 
with their typical anatomical structures and sizes. Furthermore, the boundaries of the chest wall are clearly visible 
in the dark field and the differential-phase image. 

These findings are supported by the images taken from the second mouse shown in Figure [4] While in the absorption 
image mainly the bones and the air-filled chest cavity can be seen, the other two images allow exact distinction between 
the empty chest cavity and the lung tissue. 

To evaluate the location of the lung tissue in the two mice computed tomographies (CT) of them were performed. 
This was done by rotating the object and acquiring 601 angular projections equally spaced over a full circle with the 
setup described in Section II C A reference image was taken after every 15 th projection. This leads to a total of 4.7 h 
to acquire all projection images. For image reconstruction, a filtered backprojection (FBP) was applied using a ramp 
kernel for the absorption and the dark field data and a Hilbert kernel for the phase data after performing the phase 
retrieval described in Section HTSl 

Exemplary for the data, one set of axial cross-sections of the lung region of the second mouse are displayed in 
Figure [5] (marked as horizontal lines in Figure [4]). In all three CT cross-sections, the lung tissue can clearly be located 
in the chest of the mouse. The location found with the CT reconstruction corresponds with the location visible in 
the dark field projection image and the speckled area in the differential-phase image. In contrast, the exact position 





50 100 150 200 250 300 350 400 450 




FIG. 4: Absorption (top), differential phase (middle) and dark field image (bottom) of the lung region of a second 
mouse in posterior-anterior (p. a.) orientation. In the left part of the chest cavity (right side of the images) a large 
space filled with air is visible, mainly in the dark field image. The horizontal line in each image shows the location of 

the cross-section shown in Figure [5| 



of the lung cannot be located with the absorption image alone. 



IV. DISCUSSION 

The results obtained by the investigation of in situ mouse lung tissue using grating-based X-ray phase-contrast 
imaging are very encouraging. 

First of all, three images containing different tissue properties are gained with one imaging process. This leads 
to an increase in information from which radiologists can make a more profound diagnosis. An example for the 
complementarity of information in the three images is shown in Figure [3J There, the absorption image is dominated 
by the bones of the animal and is therefore a good source for information about the skeletal structure and its diseases, 
like bone fractures. In the differential-phase image on the other hand, the bones are suppressed. Furthermore, the 
soft tissue contrast and the tissue edges are enhanced. This leads to a high organ visibility as in the differential-phase 
image of Figure [3] the lungs and the upper airways can be seen. The third image quantity, the dark field is a measure 



FIG. 5: Axial cross- sect ions of the second mouse obtained by phase-contrast computed tomography as marked in 
Figure [4] Top left: absorption image, top right: phase image, and bottom: dark field image. In this cross-section 
plane, lung tissue is only present in the right part of the mouse's chest cavity, while the rest is filled with air. 



of the granularity of structures. Because of this, the lung, which consists of many alveoli, produces a high dark field 
signal. 

Comparing the three images, it becomes obvious that the differential phase and the dark field image complement 
the absorption image and improve the identification of the shape of the lung. Especially in cases where the granularity 
of the lung tissue is affected, the dark field image can lead to a faster diagnosis and therefore treatment, as more 
information about the lung is accessible without a CT scan. It is imaginable for the identification of pulmonary 
oedema, where the alveoli are filled with water. This leads to a more homogeneous view of the lung and can, in 
combination with the higher absorption of the water visible in the absorption image, facilitate the diagnosis. Other 
examples for illnesses of the lung where the three different image quantities can possibly be exploited are lung cancer, 
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pneumonia or pneumothoraxes. For the last one, example images are shown in Figure [4] In the latter, the area which 
is occupied by the lung is just a very small part of the chest cavity. This area can clearly be identified in the dark field 
projection image, while in the absorption image the contrast between space with and without lung tissue is hardly 
visible. 

In addition to the new development of diagnostic tools for human health, this imaging approach opens a new 
avenue to possibly reduce the number of animals required for some in vivo studies. Long-term studies of effects 
of toxic substances and both soft and bony tissue may be assessed in vivo during the life of the animal. Thus, a 
significantly smaller sample size in control and verum groups may be used and long-term studies, particularly in 
academic institutions, become more feasible. 

To quantify the location of lung tissue within the mouse, a CT scan of the second mouse was performed. Its results, 
shown in Figure [5j substantiate the above results that lung tissue is, indeed, only present in the right part of the 
mouse's chest cavity in this cross-section plane. It is imaginable that an assured diagnosis can be made without a 
computed tomography scan, which could lead to a quickened diagnosis and a dose saving. 



V. SUMMARY AND CONCLUSION 

In conclusion, we presented our results on the investigation of in situ lung tissue in mice obtained with grating-based 
X-ray phase-contrast imaging. 

The three different images show different aspects of the investigated sample. While the absorption image shows 
mainly the skeletal structure, the differential-phase image enhances the tissue edges. This results in an image with 
high soft-tissue contrast and suppressed bone signal. The dark field image is sensitive to granularity in the sample. 
Therefore, the lung can easily be distinguished from the surrounding tissue. Furthermore, we demonstrated that dark 
field imaging can be used for the diagnosis of lung illnesses like a pneumothorax. For this case, we x-rayed a mouse 
with a large volume of air in the chest cavity. In contrast to the conventional absorption image, the exact lung location 
and size could easily be detected with only one projection image instead of a whole CT scan. 

In the future, we will continue our investigations on mouse tissue, not only with respect to the lung, but also other 
tissues, like the kidneys or cartilage. We will also apply this technology to investigate whether this approach can be 
applied to long-term studies in a variety of disease-associated mouse studies to provide a new diagnostic tool and, more 
importantly, to be able to reduce the number of animals needed for long-term interventions, slow healing models, and 
organ toxicity and regenerative studies. We hope to find additional cases where grating-based X-ray phase-contrast 
imaging yields an advantage compared to conventional X-ray imaging, not only to save dose, but also to increase the 
speed of an assured diagnosis. 
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